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ABSTRACT 

We report the detection of a transit egress by the ~ 3. 9- Jupiter-mass planet HD 
80606b, an object in a highly-eccentric orbit (e ^ 0.93) about its parent star of approx- 
imately solar type. The astrophysical reality of the signal of variability in HD 80606 is 
confirmed by observation with two independent telescope systems, and checks against 
several reference stars in the field. Differential photometry with respect to the nearby 
comparison star HD 80607 provides a precise light curve. Modelling of the light curve 
with a full eccentric-orbit model indicates a planet/star-radius ratio of 0.1057±0.0018, 
corresponding to a planet radius of 1.029 i?j for a solar-radius parent star; and a pre- 
cise orbital inclination of 89.285° ± 0.023°, giving a total transit duration of 12.1 ± 0.4 
hours. The planet hence joins HD 17156b in a class of highly eccentric transiting 
planets, in which HD 80606b has both the longest period and most eccentric orbit. 
The recently reported discovery of a secondary eclipse of HD 80606b by the Spitzer 
Space Observatory permits a combined analysis with the mid-time of primary transit 
in which the orbital parameters of the system can be tightly constrained. We derive a 
transit ephemeris of Ttr = HJD (2454876.344 ± 0.011) -t- (111.4277 ± 0.0032) x E. 

Key words: techniques: photometric — planets and satellites: general — planetary 
systems — occultations — stars: individual (HD 80606). 



1 INTRODUCTION 



The 'hot Jupiter' planetary companion to HD 80606, dis- 
' covered through the reflex radial-veloc ity (RV) signature 
of the parent star by iNaef et al.l (|200ir ). is notable for the 
, high eccentricity of its orbit at e ~ 0.93. The existence 
of such a system provides a laboratory for testing models 
of pla neta ry formation and ev olution (e.g., IWu fc Murray! 
(|2003h and lFord fc Rasid (|2008l ')'). particularly when the par- 
ent star is a member of a visual binary star, as here: its 
companion, H D 80607, lies at a separation of 21.1" (e.g., 
iDommanj^eFfc Nva (2Qo3 )) . corresponding to ~ 1 000 AU at 
the distance of about 60 pc l|Laughlin et al.l l|2009l ). hereafter 
L09). 

The recent observation of a secondary eclipse of 
HD80606b by L09 with the Spitzer Space Observatory gives 
insight into the effects on the planetary atmosphere of the 
highly variable insolation of radiation from the parent star 
as the planet undergoes changes in irradiation of ~ 800 fold. 

The importance of transits for characterizing the phys- 
ical and atmospheric proper t ies of exoplanets has been well 
documented (e.g. see IWinnI (|2009l )V Given the remarkable 
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nature of the orbit and irradiation of HD 80606b, a tran- 
sit observation would be especially valuable. We reiterate 
here that transits determine the orbital inclination (and 
hence true planetary mass), radius and de nsity, and thus 
constr ain internal planetary structure (see IZeng fc Seageij 
(|2008t )'). They also provide opportunities for sensitive ob- 
servation of the planetary atmospheric composition t hrough 
IR ab sorption-line studies of the transit (e.g. Tinct ti et al.l 
(120071)') ■ and emi s sion properties during secondary eclipse 
(e.g. ISwain et all (|2009l )'). Further, repeated, accurate mea- 
surement of the key parameters of the transit signal have 
the potential to probe more characteristics of the host plan- 
etary systems throug h a search f o r tim in g variations caused 
by res onant planet s (lAgol et al.l (|2005| ): Holman fc Murravl 
' 20051)). satellites llKippind (|2009l ))and Kozai migration 



Wu fc Murravl l|2003l )). Hence, the reliable detection of a 



transit of HD80606b, combined with a timing of the sec- 
ondary eclipse, would open up a very rich lode for investi- 
gation of the planet and host system. 

Following their detection of a secondary eclipse at 8fj.m 
with the IRAC instrument with Spitzer, L09 predicted the 
occurrence of a transit of HD 80606b with a-priori prob- 
ability of ~15%, according to the ephemeris (L09) Ttr = 
HJD2454653.68 + 111.4277 x iV. The predicted event of 
HJD = 2454876.5345 (2009 February 14, 0050 UT) was well 
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placed for observation from the UK, the object's dechnation 
meaning it would cross the meridian close to the zenith, and 
be available all night for continuous monitoring. Given the 
expected duration of a central transit of about 17 hours, 
and the uncertainty on the ephemeris timing, a global cam- 
paign to search for the transit ingress and egress signals was 
launched by the Transitsearch.org networkp 

In this paper, we describe a search for a predicted tran- 
sit signal by HD 80606b using the facilities of the University 
of London Observatory (ULO). Our observations with two 
independent instruments show a clear signal which we inter- 
pret to be a transit egress. 



2 OBSERVATIONS 

The University of London Observatory, operated by UCL, is 
sited in Mill Hill, NW London. Observations of known and 
candidate exoplanet transiting systems have been carried 
out for several years, and in particular have involved UCL 
undergraduate students in a campaign to monitor targets of 
specific interest. 

Several undergraduate members of the campaign were 
recruited to assist with telescope operations, and two instru- 
ments were used to monitor the field of HD 80606 through- 
out the night. 

(i) Our primary observing programme used a Celestron 
35-cm Schmidt-Cassegrain telescope on a Bisque Paramount 
ME German equatorial mount, with an SBIG STL-6303E 
CCD camera. The 22.2' x 14.8' field was centred to include 
several bright reference stars on the chip (Table [T|. An 7?- 
filter exposure time of 10 seconds was used for the entire 
series, yielding a cadence of about 20s, after readout over- 
heads. Observations of the field were acquired from 1809- 
0414 UT, and 1638 useable science frames were obtained. 

(ii) We operated a backup observing programme on a 
fork-mounted 25-cm Meade LX200 Schmidt-Cassegrain tele- 
scope, with an SBIG ST8-XME CCD camera providing a 
field of 18.6' X 12.6'. Exposure times of 30 seconds were 
used, with a typical cadence of 35 seconds. With no require- 
ment for the mount to be reversed as the target crossed the 
local meridian, these data provide an important check for 
the integrity of the Celestron photometry. A total of 761 
useable science frames were obtained through a red filter 
(Ac ~ 650 nm, AA ~ 100 nm) between 1833 and 0342 UT. 

Weather conditions were excellent throughout the 
night, and the atmospheric seeing was stable at about 2.5- 
3". The arrival of cloud shortly after 0400 UT forced all 
observations to end. 

We used standard IRAF procedures for trimming, bias 
and dark subtraction, and fiat-fielding. The bias and dark 
frames were obtained from a median stack of 10-20 frames, 
and the flat fleld from a median stack of 10 exposures of the 
evening twilight sky (Celestron) or dome (Meade). 
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Figure 1. The flux of HD 80606 relative to HD 80607. The rel- 
ative fluxes have been normalized to unity using the mean of all 
data points at HJD > 2454876.575; the Meade light curve has 
been shifted vertically by —0.023 units for clarity. The unbinned 
light curves are plotted as grey filled symbols, while J^-minute 
binned data are plotted in red. The gaps in the Celestron light 
curve at ^ .45 and ~ .51 HJD are due to a system crash and the 
reversal of the mount respectively. 



3 DIFFERENTIAL PHOTOMETRY 

We performed aperture ph otometry on a ll images using 
IRAF/DAOPHOT routines (|Stetsonl (|l987l 'l'). A tight star- 
aperture of radius 6 pixels (~ 5") and a sky annulus of inner 
radius 43" were used to avoid risk of contamination of HD 
80606 by HD 80607, and vice versa, due to their proximity. 

Fig. [T] shows the flux ratio of HD 80606 relative to the 
nearby comparison HD 80607, for both instruments. In each 
case, the difi'erential light curves have been renormalized to 
a flux ratio of unity using the mean of all data points later 
than 2454876.575 HJD. 

The light curves show a clear signal of variation, of order 
1%. Its appearance in the two independent datasets supports 
the reality of the signal: it is not an instrumental artefact. 
No attempt has been made to correct the light-curve for 
airmass-dependent effects; the observed variation is unlikely 
to be due to local, second-order, colour-dependent extinction 
effects, as the spectral types and colours of the two stars are 
very similar. 

As an additional check on the reality of the signal, we 
checked the flux ratios of both stars against the other com- 
parison stars in the Celestron field. The light curves are 
shown in Fig. [2] normalized in the same way as in Fig. [TJ 
and show a consistent variation in HD 80606 relative to all 
comparison stars. The light curves for HD 80607 are, by 
contrast, relatively flat. There does appear to be a low-level 
residual structure in the HD 80607 light curves, but we be- 
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Table 1. Properties of HD 80606 and photometric reference stars 



Star ID Sp" V' B-V^ Comparison 

HD 80606 G5 9.06 ± 0.04 0.765 ± 0.025 

HD 80607 G5 9.17 ±0.04 0.828 ± 0.029 Ref. 1 

HD 233625 G5 9.52 ± 0.04 0.62 ± 0.04 Ref. 2 

TYC 3431-0892-1 10.04 ± 0.04 1.32 ± 0.09 Ref. 3 

i g) Simbad CDS, http://simbad.u-strasbg.fr 
Hoe et al.l Il200al ). using the transformations of lMamaiek. Meyer fc LeibertI | |2002.2006|) . 



Fig. [2] lead us to choose HD 80607 as the preferred com- 
parison star for HD 80606. 
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Figure 2. The flux of HD 80606 and HD 80607 relative to com- 
parison stars Ref. 2 and Ref. 3, from the Celestron photometry. 
The relative fluxes in each case have been normalized to unity us- 
ing the mean of all data points at HJD > 2454876.575; the lower 
light curves have then been shifted vertically for clarity. 



heve these are hkely to be due to flat-field residuals caused 
by shifts of the field registration on the chip, supported by 
their apparent correlation with residuals in HD 80606 also. 

From the above, we confirm that the structure in the 
light curve seen in Fig. [T] is due to variability in HD 80606. 
There are statistical advantages to combining the fiuxes of 
all comparison stars to perform ensemble differential pho- 
tometry with higher precision in each data point from the 
increased total flux and the averaging of scintillation noise; 
however, the apparent residual systematic effects seen in 



4 LIGHT-CURVE ANALYSIS 



4.1 Fitting 



The light curves for HD 80606 illustrated in Figs. [T] 
and [2] are highly reminiscent of a transit-egress signal by 
a Jupiter-sized comp anion (e.g., see the case of HD 17156b, 
iBarbieri et al.l l|2007l ')'). In order to follow-up this hypoth- 
esis, we modelled theoretical transit light-curves using the 
orbital parameters of HD 80606b given by L0 9, a stellar 
mass of IMq (L09), and a planet mass of 3.9 Mj (|Naef et al.l 
(|200ll )'). and fitted them to the light curves presented in 
Fig [1] To account for the effe cts of high eccentricity, we 
use the dynamical model of Kipp ing (|2 008) (hcrafter, K08), 
coupled with the routines of Ma ndel fc Agoll l|200^ ') to pro- 
duce limb-darkened lightcurves. The K08 model allows fully 
for the planet's velocity to vary during the transit and for 
light-curve asymmetry. 

We combined our data into 4-min ute bins, and use d 
the AMOEBA x^-minimization routine (|Press et al.l (|l992l ')') 
to find the best-fitting light curve to the binned data. 
We assumed a solar radius for HD 80606, after L09, and 
chose quadr a tic lim b-darkening coefficients in the i?-band 
from IClaretl ||2000D for T^s = 5500 K, log 3 = 4.50 and 
[Fe/H] = 0.5, b ased on the stellar properties given by 
iNaef et all (|200ll '). 

We performed fits to three data samples: (i) Celestron 
only, (ii) Meade only, and (in) both datasets simultaneously. 
For (i) and (ii), the fits to each sample were unweighted. In- 
vestigation of the rms scatter of the residuals about each 
fit indicated that in the combined data, the Meade observa- 
tions should be weighted by 0.25 relative to the Celestron 
data. We estimated uncertainties on our fitted parameters 
by boostrap resampling of the light-curve residuals and re- 
fitting in 1000 Monte Carlo simulations. All results are given 
in Table [3 

The results of the Celestron and Meade fits are in- 
ternally consistent. Excluding the noisier data obtained at 
higher airmass at HJD < .35, the rms scatter of residu- 
als from the combined fit is 0.80 mmag (Celestron) and 
1.66 mmag (Meade). This scatter is reasonably consistent 
with the photometric uncertainties calculated from the to- 
tal fiux and scintillation noise on each star. 

The combined-fit results indicate a precise planet/star 
radius ratio of 0.1057 ± 0.0018, indicating a Jupiter-sized 
companion for a stellar radius of 1 _Rq. The inclination is 
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Figure 3. The Jf-minute binned data for the Celestron (filled 
symbols ) and Meade ( open ) differential light curve of HD 80606 
with respect to HD 80607, with the best-fitting eccentric- orbit light 
curve drawn through the points. The residuals from the best-fit 
curve are shown, at a relative flux of +0.975 units. 



tightly constrained at 89.285±0.023 degrees, yielding a total 
transit duration, Ti^4, of about 12 hours. 

4.2 Consistency of the signal with the predicted 
transit of HD 80606b 

Since the transit-like signal occurs in both data sets, we 
conclude it is of astrophysical origin. HD 80606 itself is not 
expected to be variable: indeed, the chromospheric activity 
of the star as measured b y its Can H and K core emission 
jSaffe. Gomez fc Chaverol |2005)) is indicated to be quiet. 

We can show that the mid-transit time of our best-fit 
model is consistent with the ephemeris of L09. Using the 
model of K08, we calculate the difference between the pri- 
mary and secondary eclipse mid-points. The orbit param- 
eters and uncertainties of L09 then give a predicted mid- 
transit of Tmid = 2454876.44 ± 0.24 HJD, compared to 
our observation of 2454876.344 ± 0.011 HJD, ~ 2.5 hours 
early but well within the predicted uncertainty. The ob- 
served egress duration and depth constrains the second con- 
tact to occur at ~ 1650 UT on 2009 February 13th, con- 
sistent with the absence of an observed ingress. Further, 
we can use the K08 model with our best fit parameters to 
predict the secondary-eclipse duration. Our model predicts 
Si,4 — 1.93 hours and 82,3 — 1.56 hours, completely consis- 
tent with the results of L09 who found si,4 — 1.92 ± 0.10 
and S2,3 = 1.68 ± 0.24 hours. 

We therefore conclude that there is a very high proba- 
bility that the observed signal is indeed a planetary transit 
of HD80606b. 
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Figure 4. The vertical lines represent the best-fit, one-sigma and 
two-sigma confidence limits on e and vo. The black ellipse repre- 
sents the previously derived one-sigma constraints from LOO. The 
two limits combined provide a highly precise measurement of e. 



4.3 Eccentricity constraints 

Based on the timing between the secondary and primary 
transit, it is possible to derive precise constraints on the 
orbital eccentricity. L09 cited that the primary transit occurs 
only ~ 6 days after the secondary eclipse, short in relation 
to the 111-day period due to the extremely high eccentricity. 
The eccentricity has been derived from RV measurements, 
and the time delay between primary and secondary eclipse 
can be used to constr ain it very t i ghtly. Typically, previous 
authors, for example IWinn et al.l (|2005l ). have used a first- 
order approximation of e for such a calculation to derive 
e cos t*7, but in this case we are very far from a 'low' value of 
e. We have therefore used the K08 model to vary e and w 
such that it reproduces the observed delay of 5.897 ± 0.017 
days and we plot the allowed parameter values in Fig. 4. 

Our constraints on e and 117 can be seen to intersect 
the 1-sigma parameter space of L09 derived from indepen- 
dent RV measurements. Combining both constraints, we 
derive highly precise values of e = 0.9336 ± 0.0002 and 
■cu = 300.4977 ± 0.0045, or a precision in the eccentricity 
of 2 parts in 10000. 



5 DISCUSSION AND CONCLUSIONS 

We have detected a transit-egress signature towards HD 
80606 in two sets of photometry taken simultaneously on 
the night of February 13/14 2009 from the ULO, London, 
UK. The detection of the signal in two independent sets of 
photometry strongly supports an astrophysical origin for the 
signal. Further, the moment of mid-transit is consistent with 
the ephemeris of L09 as well the egress duration being con- 
sistent with both the absence of observed ingress and the 
measured secondary-eclipse durations of L09. We conclude 
that the signal is highly likely to be that of the transit of 
HD 80606b. 

The planetary radius of 1.03 R,j, for a stellar radius of 
IRq implies a planet density of 4440 ± 240 kg m~^, surpass- 
ing the density of the other known highly eccentric transit- 
ing planet, HD 17156b, and making HD 80606b an object 
of great interest to the field of planetary formation. 
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Table 2. Fitted planetary parameters 





C-14 


Meade 


Combined 


Number of binned data 


139 


134 


273 


Central transit depth (%) 


1.048 


1.116 


1.059 


Ratio-of-radii, k 


0.1043 ±0.0017 


0.1100 ±0.0034 


0.1057 ±0.0018 


Egress duration, T3,4 (hours) 


2.84 ±0.22 


3.3 ±0.8 


2.98 ±0.26 


Total duration, Ti^4 (hours) 


12.3 ±0.4 


11.9 ±0.6 


12.1 ±0.4 


Mid-transit time, T^\^ 


0.339 ±0.011 


0.349 ± 0.018 


0.344 ±0.011 


(HJD - 2454876) 








Orbital inclination, « (°) 


89.296 ±0.032 


89.27 ±0.05 


89.285 ± 0.023 


Planetary radius, Rp (Rj) 


1.015 ±0.016 


1.070 ±0.030 


1.029 ±0.017 



Hence, HD 80606b has the longest period and most ec- 
centric orbit amongst the known transiting planets, making 
it an extremely unusual and interesting object for further 
study. The large change in stellar irradiation between pri- 
mary and secondary eclipse offers an excellent opportunity 
for follow-up atmospheric studies. 

We note that very long transit duration allows for ex- 
tremely precise measurements of mid-transit time, depth 
and inclination. Therefore, any periodic or secular modifi- 
cation to these parameters due to a perturbing body would 
be easier to detect relative to other systems. The planet's 
orbit is much wider than other transiting planets, and tak- 
ing into account a significant light-time travel effect of 2.75 
minutes between secondary and primary eclipse, we have 
demonstrated that the eccentricity of this system can be de- 
termined to a level of 2 parts in 10000. We derive a transit 
ephemeris of Ttr = HJD (2454876.344±0.011) + (111.4277± 
0.0032) X E. 

Finally, the work reported here continues to highlight 
the utility of co-ordinated global campaigns which promote 
the use of relatively modest telescope apertures to observe 
transits and help characterize exoplanet systems; indeed, 
such campaigns have been actively encouraged through the 
oklo.org website of G. Laughlin. The full characterization of 
the transit signal of HD 80606b will almost certainly depend 
on future co-ordination of observations made from multiple 
sites. We also emphasize the opportunities for useful work to 
be carried out from observing sites in locations — even Lon- 
don — which are not usually considered to be suitable for 
modern observational astrophysics. While there are bright 
stars hosting known or candidate transiting planets await- 
ing detection or characterization, this will remain true. 

[We note that at the time this paper was submitted, two 
independent reports of a detection of this same transit event 
were announced — see Moutou et al. ( arXiv:0902.4457 J and 
Garcia-Melendo & McCuUough ( arXn^QOOWg]) .] 
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